Astronomy & Astrophysics manuscript no. draft021737 


© ESO 2009 


February 27, 2009 





The mass profile of early-type galaxies in overdense environments: 
the case of the double source plane gravitational lens 

SL2SJ021 76-051 3* 

H. R. Gavazzi2'3 4, M. Limousin^ '^ ^ R. Cabanac^ R J. Marshall^, B. Fort^-^ T. Treu^, R. Pello^ E. Jullo^-'', 

J.-R Kneib^''°, J.-F. Sygnet^-^ 

' Physics Department & Shanghai Key Lab for Astrophysics, Shanghai Normal University, 100 Guilin Road, Shanghai 200234, 
China 

e-mail: tuhongdshnu . edu . cn 
^ CNRS, UMR7095, Institut d'Astrophysique de Paris, 98bis Bd Arago, F-75014, Paris, France 
3 UPMC Univ. Paris 6, UMR7095, Institut d'Astrophysique de Paris, 98bis Bd Arago, F-75014, Paris, France 

* Department of Physics, University of California, Santa Barbara, CA 93 106, USA 

' Laboratoire d'Astrophysique de Toulouse-Tarbes, Universite de Toulouse, CNRS, 57 avenue d'Azereix, F-65000 Tarbes, France 

* OAMP, Laboratoire d'Astrophysique de Marseille - UMR 61 10 - Traverse du siphon, F- 130 12 Marseille, France 
^ Dark Cosmology Centre, Niels Bohr Institute, University of Copenhagen, 2100 Copenhagen, Denmark 

^ Laboratoire d'Astrophysique de Toulouse-Tarbes, CNRS-UMR 5572, and Univ. Paul Sabatier, F-31400 Toulouse, France 
' European Southern Observatory, Alonso de Cordoba, Santiago, Chile 
Department of Astronomy, California Institute of Technology, 105-24, Pasadena, CA91125, USA 



Preprint online version: February 27, 2009 

ABSTRACT 

Context. The Strong Lensing Legacy Survey (SL2S) provides a sample of strong lensing events associated with massive distant 
galaxies, some of which lie in the outskirts of galaxy groups and clusters. 

Aims. We investigate the internal structure of early-type galaxies in overdense environments, where tidal forces are expected to alter 
dark matter halos of infalling galaxies. 

Methods. SL2SJ02176-0513 is a remarkable lens for the presence of two multiply-imaged systems at different redshifts lensed by a 
foreground massive galaxy at zicna = 0.656: a bright cusp arc at Zmc = 1-847 and an additional double-image system at an estimated 
redshift of zm ~ 2.9 based on photometry and lensing geometry. The system is located about 400 kpc away from the center of a 
massive group of galaxies. Mass estimates for the group are available from X-ray observations and satellite kinematics. Multicolor 
photometry provides an estimate of the stellar mass of the main lens galaxy. The lensing galaxy is modeled with two components 
(stars and dark matter), and we include the perturbing effect of the group environment, and all available constraints. 
Results. We find that classic lensing degeneracies, e.g. between external convergence and mass density slope, are significantly 
reduced with respect to standard systems and infer tight constraints on the mass density profile: (i) the dark matter content of the main 
lens galaxy is in line with that of typical galaxies /dn,(< ^c) = 0-41^[5 q6; (ii) ths required mass associated with the dark matter halo of 
the nearby group is consistent with X-ray and weak-lensing estimates (cr^p = 550*™)^', (iii) accounting for the group contribution in 
the form of an external convergence, the slope of the mass density profile of the main lens galaxy alone is found to be a = -1.03;^q j^, 
consistent with the isothermal (a = -1) slope. 

Conclusions. We demonstrate that multiple source plane systems together with good ancillary dataset can be used to disentangle 
local and environmental effects. 
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1. Introduction 

In the last thirty years, gravitational lensing has become a ma- 
jor astrophysical tool because of its unique ability to probe the 
densest regions of the universe independent of their emitted radi- 



* Based on observations obtained with MegaPrime/MegaCam, a joint 
project of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii 
Telescope (CFHT) which is operated by the National Research Council 
(NRC) of Canada, the Institut National des Science de I'Univers of the 
Centre National de la Recherche Scientifique (CNRS) of France, and 
the University of Hawaii. This work is based in part on data products 
produced at TERAPIX and the Canadian Astronomy Data Centre as 
part of the Canada-France-Hawaii Telescope Legacy Survey, a collabo- 
rative project of NRC and CNRS. 



ation. Large imaging or spectroscopic surveys are becoming key 
players, in that context, by bringing numerous lensing events of 
all scales and spanning an increasingly broader range of masses 
and redshifts (e.g. Kochanek et al. 1999; Myers et al. 2003; 
Inada et al. 2008; Bolton et al. 2006, 2008a; Faure et al. 2008a; 
Marshall et al. 2008). 

Among the imaging searches for lenses, the CFHT Legacy 
Survey' stands out by combining depth in five photometric 
bands, excellent image quality (~ 0'.'7 - 0'.'8 median seeing) 
and an area of 170 deg^ when completed. The Strong Lensing 
Legacy Survey (SL2S)^ aims at automatically detecting grav- 



' http://www.cfht. Hawaii .edu/Science/CFHLS/ 
- http://www-sl2s.iap.fr/ 
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itational lensing events in these data (Cabanac et al. 2007; 
Limousin et al. 2008), be they either large arcs in groups and 
clusters, or smaller, galaxy-scale lensing events (sometimes also 
referred to as "Einstein rings" to remind us that the lensed back- 
ground source is extended, and therefore does, at some surface 
brightness level, produce a complete ring around the deflector). 
A detailed description of the ringf inder detection method will 
be given in a forthcoming paper (Gavazzi et al. 2009, in prep). 

One of the scientific questions that can be addressed in a di- 
rect and powerful manner by gravitational lensing is that of the 
dark matter halo profiles of early-type galaxies, the most com- 
mon type of strong lenses. For example, gravitational lensing 
has been used to demonstrate the presence of dark matter ha- 
los (e.g. Kochanek 1994; Treu & Koopmans 2004), and to test 
the vaUdity of the universal dark matter halos predicted by nu- 
merical simulations (Navarro et al. 1997a; Gavazzi et al. 2007). 
When galaxies fall into an overdense environment, i.e. a cluster 
or a group, their dark matter halos are expected to be stripped 
by tidal interactions (e.g. Dobke et al. 2007). Tidal stripping 
typically does not reach deep enough into the infalling galaxy 
to afl'ect its luminous structure (e.g Treu et al. 2003) until the 
galaxy reaches the densest regions, and therefore the effects are 
hardly detectable with classical astronomical probes. However, 
stripping of the outer halo begins much earlier during the in- 
fall and can be detected by gravitational lensing, for example 
by weak lensing analysis of cluster galaxies (Natarajan & Kneib 
1997; Geiger & Schneider 1998; Natarajan et al. 2002a, 2007). 
Analyses of selected galaxy samples coupling lensing and dy- 
namics have also suggested that halos of satelHte galaxies may 
be truncated and efl'ectively steeper than those of isolated galax- 
ies (e.g. Treu & Koopmans 2002; Auger 2008; Treu et al. 2008). 
In turn, galaxy-scale strong lenses in the field of clusters, can be 
used to constrain the cluster mass distribution near the high mag- 
nification regions of the cluster (Tu et al. 2008; Limousin et al. 
2008) and that of nearby structures (e.g. Kochanek & Blandford 
1991; King 2007; Tu et al. 2008). 

Beyond the general problem of tidal truncation of dark mat- 
ter halos, the envirormient of lens galaxies has been a topic of in- 
tense debate in the lensing literature for a variety of reasons. One 
important question is the role of the environment in the measure- 
ment of Hubble's constant from lens time delays (e.g. Kochanek 
2002; Kochanek & Schechter 2004). In addition, the environ- 
ment of lens galaxies must be understood if one wants to gen- 
eralize the findings to the overall population of similar galaxies. 
It is known that multiply-imaged QSOs are often found behind 
overdense regions which have noticeable efl'ects on the lensing 
configurations (e.g. Keeton et al. 1997; Fassnacht & Lubin 2002; 
Auger et al. 2007; Momcheva et al. 2006; Oguri 2006; WiUiams 
et al. 2008). This preference for dense environments is also seen 
in the low-redshift SLACS sample (Treu et al. 2008). However, 
these authors showed that this preference is no more pronounced 
than for any massive elliptical galaxy, whether or not it is ac- 
tually a lens. Within the errors, the preference for high-density 
environments is just a consequence of the clustering of massive 
galaxies. A somewhat similar trend is observed in the COSMOS 
survey (Faure et al. 2008b) although the use of weak lensing 
convergence map as a proxy for the environment makes the in- 
terpretation a bit more difllcult. This preference for clustering 
is nicely illustrated by the results of Fassnacht et al. (2006) and 
Newton et al. (2008) which can be summarised as: a good place 
to look for new lenses is the vicinity of known lenses. 

In this paper, we present the analysis of an extraordinary sys- 
tem SL2SJ02 176-05 13 discovered by the CFHT Legacy Survey, 
and use it as a case study to investigate the eff"ect of environ- 



Table 1. Summary of the most relevant observational quantities 
for SL2SJ02176-0513. The eUipticity is defined as e = (a^ - 
lP')l(d^ + b^). Masses are given with their photometric errors. 
References: " Geach et al. (2007). Note that errors on photometry 
are dominated by systematic calibration uncertainties. 



Lens galaxy centre (J2000) 


RA 


02:17:37.136 




DEC 


-05:13:29.41 


Lens galaxy redshift 


Zlens 


0.6459 ± 0.0003 


Lens mags {R < 1") 


U* 


23.90 ± 0.20 




g' 


22.89 ± 0.07 




r' 


21.74 + 0.05 




i' 


20.44 ± 0.04 




z' 


20.12 ±0.03 


HST (Sersic model mags) 


F606W 


20.90 ± 0.05 


F814W 


19.71+0.01 


Effective radius (F814W) 


R. 


a'90±a'io 






6.20 ± 0.7 /i-'kpc 


Sersic model index (F814W) 


n 


5.20 


Rest frame V band luminosity 


Lv 




light position angle (F814W) 


PA, 


78 ± 6 deg 


light ellipticity(F814W) 


e. 


0.14 + 0.06 


"Photometric" stellar mass 


M, 


2.1-4.1 X 10" Mo 


Tangential arc 






. . . redshift 


Zarc 


1.8470 + 0.0003 




F606Waic 


22.04 + 0.16 




F814Waic 


22.19 + 0.13 


Double system (Rl) 






. . . photom. redshift 


Zdbl 






F606Wri 


24.11+0.06 




F814Wri 


23.91 + 0.06 


Einstein radius of the main arc 


Rare (") 


1.4 


Group redshift" 


Zgrp 


0.648 + 0.001 


Group velocity dispersion" 


spec 
CTgrp 


770+ 170 km s"' 


X-ray temperature" 


kT 


2.0:i:«keV 


cr from cr-Tx relation" 




520+ 120 km S-' 



ment on the dark matter halos of galaxies falling into over- 
dense environment. The system is unusual for several reasons. 
First, it is producing multiple images of two sources at different 
redshifts, similarly to the double Einstein Ring discovered by 
Gavazzi et al. (2008). Second, the system is located near a well- 
constrained group environment, and we were able to gather a 
large spectrophotometric set of information for both the lensing 
galaxy and its environment. Using this wealth of information, 
we are able to break some of the classic degeneracies of lensing 
studies. This allows us to infer the relative distribution of Ught 
and dark matter in the lens galaxy out to 1 .5 effective radii of the 
lens galaxy, and probe the effect of the nearby group potential. 

The paper is organised as follows. After presenting the data 
set available for SL2SJ02 176-05 13 in Sect. 2, we measure the 
optical properties of the lens galaxy (including stellar mass) in 
Sect. 3, and then present our lens modelHng in Sect. 4. In Sect. 6, 
we discuss the main conclusions from the modelling and future 
possible improvements. 

Unless otherwise stated we assume a concordance cosmol- 
ogy with Hq = IQhjo kms"' Mpc"', = 0.3 and Qa = 0.7. 
At the redshift of the lens SL2SJ02 176-05 13 zie„s = 0.6459, 
1" ~ 6.9 kpc. All magnitudes are expressed in the AB system. 



2. Observational data 

In this section we present the ensemble of available data. Table 1 
summarises the main quantities of interest for this work. 
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Fig. 1. Wide field CFHT image of SL2SJ02 176-05 13 (irg "inverted" composite image). The big 10" radius circle shows the location 
of the nearby galaxy group studied by Geach et al. (2007). 



2. 1 . SL2SJ021 76-0513 in the CFHTLS imaging data 

The SL2S detection strategy of lensing events was introduced by 
Cabanac et al. (2007). Details of the ringfinder detection al- 
gorithm as well as the full sample of gravitational Einstein rings 
will be presented in a forthcoming paper (Gavazzi et al. 2009, in 
prep.). As part of this effort the gravitational lens SL2SJ02176- 
0513 was detected as a promising galaxy-scale lens candidate. 
Subsequent high-resolution follow up imaging and spectroscopy 
were since obtained. SL2SJ02 176-05 13 belongs to the Wl field 
in the wide part of the CFHTLS survey. Currently, sub-arcsecond 
seeing imaging is available in the Megacam bands u*, g', r', 
i' and z'. A large scale colour image of SL2SJ02 176-05 13 is 
shown in Fig. 1. The coordinates of the elliptical lens on the 
CFHTLS image is ay20oo=02:17:37.136, 5y20oo=-05;13;29.41. 



ies (see Fig. 5 in Geach et al. (2007)) shows that the measured 
group velocity dispersion is somewhat higher than one would 
expect from its X-ray luminosity: the Lx - cr relation predicts 
cr^p - 520 ± 120 km s"'. We can notice on Fig. 6 of Geach 
et al. (2007) that the positive and negative radial velocities of the 
group members do not seem to be distributed randomly across 
the group, but are rather distributed in a symmetric way on each 
side of the X-ray group centre. This may be a hint that the group 
is not relaxed, and hence that the measured spectroscopic ve- 
locity dispersion might best be taken as an upper limit on the 
depth of the potential well. We note that Geach et al. (2007) also 
identified the lensing nature of the system and measured both 
the lens and doubly-imaged background source redshifts, find- 
ing ziens - 0.646 and z^ic - 1.847. 



2.2. Radio and X-ray properties 

Simpson et al. (2006) found this object to be a radio source 
(VLA J2000 02: 17:37.21 -05: 13:27.96 ; SXD J2000 02: 17:37.21 
-05:13:28.0). Its environment was investigated by Geach et al. 
(2007). The lens galaxy inhabits a massive galaxy group (or 
equivalently a low mass galaxy cluster) at Zgrp = 0.648 + 0.001, 
with an optical spectroscopic velocity dispersion CTg," - 770 + 
170 km s"' . They also quote a group X-ray temperature of kTx - 
1 .97;^p"gkeV, and a group X-ray luminosity of Lx, [OJ.iokeV] - 
1.79!°-g X lO'^^erg/s. Geach et al. (2007) give the XMM map 
that allows us to set the centroid of the group extended X-ray 
emission to be Q'y2()00=02: 17:33.534, 5y2ooo=-05:13:14.93, with 
an accuracy of +10": the group centre is therefore 56 + 10" 
from SL2S J02 1 76-05 1 3 . The Lx-cr relation for groups of galax- 



2.3. HST observations 



SL2SJ02 176-05 13 was observed with the Hubble Space 
Telescope Advanced Camera for Surveys (HST/ACS) in the 
F814W and F606W bands, as part of our SL2S follow-up snap- 
shot program (SNAP10876, PI Kneib). A standard reduction was 
done with multidrizzle and we used LACosmic for cosmic 
ray removal. This method happened to be more efficient in the 
F814W band image which is made of 2 x 400s exposures, but 
gave satisfying results in the single 1 x 400s F606W exposure. 
Fig. 2 shows an image of SL2SJ02176-0513 in this latter band, 
with lensed features labelled (see Sect. 4.1). 
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Fig. 2. HST image (F606W band) of SL2SJ02 176-05 13. The 
lens redshift is ziens = 0.6459. The main conjugate points 
Q1,Q2,Q3,Q4 forming the cusp arc (at redshift Zm-c = 1.847) are 
shown as well as the additional double-image system R1-R2. 

2.4. Keck I spectroscopy 

We obtained a spectrum of SL2SJ02 176-05 13 with LRIS spec- 
trograph on the Keck I telescope during the night of December 
23rd 2006. Observations were made in long slit mode (r.'5 wide) 
and oriented East- West. The spectrum therefore presented two 
parallel tracks that allowed us to measure at the same time the 
lens galaxy redshift ziens = 0.6459 + 0.0003 (based on the typi- 
cal absorption features like Call H & K, G band, Mgb) and the 
lensed arc redshift Zm-c = 1 .847, using the strong emission lines 
Lya, CIV, Hell, OIII], and CIII]. Most of them are typical of a 
narrow line type II AGN (e.g. Steidel et al. 2002). This suggests 
that the radio source and associated active nucleus, identified by 
Simpson et al. (2006) and Geach et al. (2007), might actually 
not be at the group redshift Zgip ^ ziens but rather at redshift 
1.847. High resolution VLBI imaging would be able to demon- 
strate whether the radio emission is consistent with being lensed 
by the foreground galaxy. 

3. Lens galaxy properties 

The HST imaging allows us to make precise shape measure- 
ments of the lens galaxy optical surface brightness. Using the 
galf it facility (Peng et al. 2002), we found that the lens galaxy 
is well fit by a Sersic surface mass density profile (assuming 
constant stellar mass-to-light ratio) of the form 

m) = 2, exp [-b„ {(R/R,)"" - l)] , (1) 

with Sersic index n ^ 5.2 (corresponding to b„ ^ 10.07) and 
effective radius R^.. We thus measured total model magnitudes 
F606W = 20.90 + 0.05 and F814W = 19.71 + 0.01 and effec- 
tive radii R^ = 0'.'90 + a.'lO ^ 6.20 ± 0.7 h'^ kpc in the F814W 
band. In addition, galfit yielded the axis ratio = b/a - 
0.87 + 0.05, such that the ellipticity e, = 0.14 ± 0.06 where 
e = (1 - q^)/(l + q^), and the position angle PA. = 78 ± 6° 
(North to East, ccw)"*. The measured flux in F814W band was 
converted into an absolute rest frame V band luminosity of 
Ly - (1 .36 ± 0.08) X 10" /lyQ Lq. This includes corrections from 
Galactic extinction (Schlegel et al. 1998), the small filter shift- 
ing correction to translate z = 0.6459 F814W into a rest frame 

^ Note that systematic uncertainties dominate the galfit run. We 
use the difference between estimates in the F814W and F606W bands 
as an error estimate for both q, and PA, 



V band, in nearly the same way as (Gavazzi et al. 2007) except 
that, here we do not apply any passive evolutionary correction to 
the reference redshift zero epoch. To do so, one could use the re- 
lation constrained by Fundamental Plane evolution studies (Treu 
et al. 2001,2006): 



CFHT photometry within a 1" radius aperture around the 
lens centre allows us to estimate the corresponding cylindri- 
cal stellar mass. Note that we applied small ~ 4% corrections 
to compensate aperture flux losses due to ground-based seeing 
(ranging from I'.'l in u* to 0'.'6 in i'). CFHT magnitudes were 
also corrected for Galactic extinction. Using the publicly avail- 
able code hyperz (Bolzonella et al. 2000) with a Chabrier Initial 
Mass Function (IMF) (Chabrier 2005) and an intrinsic absorp- 
tion coefficient Ay - typical for nearby early-type galaxies, 
we obtain a stellar mass ofM,(< 1") = 1.4+"-2 x 1O"M0. Using 
a Salpeter IMF and/or a non-zero A v up to 0.3 (unlikely for such 
a massive early-type galaxy), the inferred mass reaches Mt{< 
1") ^ 1.7 x1O"M0. It results that, although dominated by stellar 
evolution systematics, the stellar mass of the lens galaxy within a 
one-arcsec radius is likely in the range 1.2-2.3xlO"M0. Using 
the Sersic profile measured on HST data, we can extrapolate this 
aperture mass measurement and estimate a total "photometric " 
stellar mass M» = 2.1 -4.1 x 1O"M0. This value corresponds to 
a V-band steflar Mass-to-light ratio MJLy - 1.5-3.0 (M/L)q at 
Ziens = 0.646. Accounting for luminosity evolution as described 
above (Eq. 2), we get an equivalent zero-redshift M./Ly = 
2.8 - 5.5 {M/L)q, which is consistent with, although slightly 
lower than, local estimates (Gerhard et al. 2001; Trujillo et al. 
2004), but is in good agreement with the strongH-weak gravita- 
tional lensing constraints of Gavazzi et al. (2007), who found an 
average M,/Ly - 3.84 ±0.40 in SLACS lens galaxies, once they 
were evolved to redshift zero using Eq. (2). 

4. Strong Lensing Analysis 

In this section, we present our strong lensing analysis. It was 
performed using the latest version of the lenstool code (Julio 
et al. 2007), which uses a Markov Chain Monte Carlo sampler 
to characterise the posterior PDF for the model parameters. 

4.1. Identification of multiple images 

lenstool relies on the identification of conjugate multiple 
images of the same background source. For a source imaged 
times, since the unlensed source position is unknown, the loca- 
tion of each of these points in the image plane results in 2{N - 1) 
independent constraints. When a source is structured in such a 
way that multiple sub-components (surface brightness peaks) 
can be identified, one can significantly increase the amount of 
information provided by the lensing configuration. 

For the HST/ACS images we are considering here, we find 
the typical astrometric uncertainties on the positions of identi- 
fied conjugate image features to be 0'.'03. These constraints are 
summarised in Table 2, while the subsections below detail the 
way in which multiple images are identified. 

The HST/ACS image in Fig. 2 clearly reveals that 
SL2SJ02 176-05 13 is causing the formation of a tangential arc 
(Q1-Q2-Q3, in a typical cusp configuration) along with its 
small and faint counter-image (Q4) on the opposite side of the 
lens, about twice as close to the lens centre. The arc is bright 
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Table 2. Positions of multiply-imaged, conjugate image features 
used in the lenstool analysis. Units are arcseconds West and 
North relative to the lens galaxy centre. Uncertainties are set to a 
typical value of (X'OS. Column "# const" gives the contribution to 
the total amount of constraints supplied to the lenstool model. 



Image 1 Image 2 Image 3 Image 4 # const 

Q -1.06,1.03 -1.29,-0.53 -1.06,-0.98 0.88,0.08 6 

Q' -1.33,0.41 -1.34,-0.25 -0.90,-1.18 0.88,0.08 6 

R -1.28,2.31 0.39,-0.33 2 

R' -1.50,1.89 0.45,-0.25 2 

total 16 



are predicted by our lens model to form well underneath the el- 
liptical galaxy, thus explaining why we do not observe them. 
This non-detection prevents us from accurately fixing the loca- 
tion of the radial critical line, but the relative distances between 
Rl-R'l and R2-R'2 still allow us to extend the range of radial 
constraints. 



(F606W,,c = 22.04 ± 0.16 and F814Warc = 22.19 ± 0.13). A 
rough estimate of its length-to-width ratio r yields r ^ 8, which 
means that the unlensed source magnitudes are ~ 2.3 mag dim- 
mer The arc has a well-extended surface brightness distribution 
resulting from the merging of 3 conjugate multiple images sep- 
arated by the critical line. The source thus has to lie across the 
caustic with parts of the source pointing outside the tangential 
cusp. Those parts of the source will only be imaged 3 times, 
whereas the others will be imaged 5 times^. This configuration 
makes the identification of conjugate images of substructures 
along the arc quite difficult ex nihilo. 

We approach this problem in two stages. In the first stage, 
we perform a modelling using the most readily identified (and 
therefore most robust) conjugate points in the arc and its counter- 
image only. This provides a satisfactory model which we then 
use for testing further the more difficult conjugate images along 
the arc, finally reaching the astrometric precision of 0'.'03 in the 
predicted image points. 

An interesting feature of SL2SJ02 176-05 13 is that a highly- 
elongated arclet shape is present at about 2'.'5 North-East of the 
lens galaxy (labelled Rl in Fig. 2). A photometric redshift of this 
system Zdbi - 2.90^q 24 was obtained from the CFHTLS pho- 
tometry with hyperz. The lens modelling of the tangential sys- 
tem Q1-Q2-Q3-Q4 predicts that Rl, at that redshift, should be 
multiply-imaged. Indeed, we do observe a demagnified counter- 
image R2 (clearly in the F606W image, and at slightly lower 
signal-to-noise in the F814W image). R2 is unfortunately too 
faint and buried in the light of the lens galaxy for an accurate 
confirmation of its photometric redshift. However, in the HST 
imaging we were able to measure photometry for R2 and Rl, 
leading to colour indices (F814W-F606W)ri = 0.19+0.08 and 
(F814W- F606W)r2 = -0.15 ± 0.85 that are at least consistent 
with one another, although this constraint is relatively poor due 
to the extreme faintness of R2 (F606W = 27.15 + 0.54; Rl is 
much brighter, F606W = 24.1 1 ± 0.06). 

We can also estimate the "geometric redshift" of the R1-R2 
source by including it in the lens modelling. We identify two sec- 
ondary conjugate brightness peaks R'l-R'2, and use these with 
R1-R2 to provide 4 new astrometric constraints. We infer the 
value of Zdbl to be very close to the photometric value. 

Throughout the rest of this work, we fix the redshift Zdbi 
of the double-image system to its photometric redshift estimate 
Zdbl = 2.9, noting that the uncertainty on this value is small 
enough for the corresponding error in the Dx^jD^ distance ratio 
to be negligible. In addition one might expect a third pair of im- 
ages R3-R'3 of opposite parity to R2-R'2 with respect to the ra- 
dial critical line - if it exists. These demagnified counter-images 

In practice, this odd number of conjugate images reduces to 2 and 
4 multiple images respectively because one of the images (the central 
one) is highly demagnified and almost never observed 




Fig. 3. Illustration of the accuracy of the best fit model. Small 
blue circles are input data while red crosses are model-predicted 
image positions after optimisation. Critical curves are shown as 
red solid lines. For reference, the distance of the galaxy centre to 
the main arc is about 1'.'4. 



4.2. Model assumptions 

We now describe our mass model for SL2SJ02 176-05 13 and its 
group environment. The lens itself is modeled with two com- 
ponents describing the stars in the lens galaxy, and the dark 
matter (DM) halo in which they sit. For the surface density of 
the former, we assume a Sersic profile for which most model 
parameters (effective radius Rg, index n, ellipticity and orienta- 
tion) are fixed at the values determined from the galfit model 
of the HST images. Only the total stellar mass M» is taken as 
a free parameter. Given the total rest-frame V band luminosity 
Lv - (1.36 + 0.08) X 10" /i^Q Lq, one can readily translate this 
parameter into a stellar mass-to-light ratio. 

The dark matter halo surrounding the lens galaxy is modeled 
as a Pseudo Isothermal Mass Distribution (PIEMD; Kassiola & 
Kovner 1993; Limousin et al. 2005). We also considered a NFW 
profile (Navarro et al. 1997b) as suggested by N-body simula- 
tions, and found that none of our results change significantly ei- 
ther in terms of goodness-of-fit or regarding the main qualitative 
trends described in the results section below. We fix the centre 
of the halo to that of the galaxy light distribution, building on 
well-established results from previous galaxy-scale strong lens- 
ing studies (e.g. Yoo et al. 2006). (Moreover, when the centroid 
was freed in our modelling, we found that the centre of the DM 
halo has to be < O'.'l from the centre of light.) The orientation 
PAdm, ellipticity edm, and velocity dispersion parameter cTdm of 
the DM halo are let free with broad uniform priors < edm ^ 0.6, 
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< PA±n < 180°, and <Tdm < 450 km s''. We also allow the 
possibility of a finite core radius in the range < r^^^ < I'.'S. 
Due to the relatively small projected distance to the cluster cen- 
tre 400 kpc), the galaxy may have a large part of its DM halo 
stripped. Lacking observational constraints, we set the trunca- 
tion radius rcut to 60 kpc, an average value motivated by simula- 
tion and galaxy-galaxy lensing results in similarly high density, 
cluster environments (Natarajan et al. 2002b,c; Limousin et al. 
2007a). 

Finally, the host halo of the galaxy group to which 
SL2SJ02 176-05 13 belongs will have some impact on the lens- 
ing properties of the main lens galaxy. To first order, the per- 
turbing effect of the group can be characterised by an external 
shear, oriented tangentially relative to the group centre. A sub- 
stantial amount of extra convergence should also be supplied by 
this nearby group DM halo. Since we know where the group is 
(from the X-ray observations. Sect. 2), we thus decide to model 
the group by a Singular Isothermal Sphere (SIS) with veloc- 
ity dispersion cTgip. This model not only predicts the external 
shear and its associated convergence, it also naturally accounts 
for any higher order deflections at the lens position. Typically, 
if the group velocity dispersion suggested by the kinematics of 
group member galaxies or the intra-cluster gas temperature is 
correct (cTgjp ~ 630 km s~'), one would expect the group to pro- 
duce external shear and convergence of order lyextl ~ fext - 0.05 
at the location of the lens. We thus include o-gyp as the last free 
parameter of the model, assigning a uniform prior < o-g^ < 
900 km s~' in order to be able to compare the group mass esti- 
mate with the independent X-ray and spectroscopic values. 

Altogether we have a total of 6 free parameters in our lens 
model. However, instead of further considering the raw model 
parameters above, we introduce the following secondary trans- 
formed quantities, more relevant for the discussion below: the 
group velocity dispersion o-g^p, the total stellar mass of the lens 
galaxy M,, the dark matter fraction enclosed within the cylinder 
defined by the effective radius /dm(< ^e). the difference in ori- 
entations of the dark matter halo and the stellar component APA 
(for which the errors in the light orientation PA^ = 78 ± 6 are 
added in quadrature to the PA^^ estimate) and the ellipticity of 
the dark matter halo e^. 



5. Results 

The lens modelling yields good fits to the positions of all the 
multiple images. The recovered model parameters, along with 
their 68% confidence level uncertainties and goodness-of-fit 
- quantified both in terms of per degree of freedom and 
Bayesian log-Evidence - are listed in Table 3. In Fig. 3, we 
show the excellent matching between the observed positions of 
the multiple images (from Table 2) and the ones predicted by the 
best-fit model. We stress here that the combination of medians 
of marginalised one-dimensional posterior distribution for each 
parameter (as listed in Table 3) does not necessarily correspond 
to the mode of the overall posterior probability distribution that 
we nominate "best fit model" and approximate with the highest 
Hkelihood MCMC sample. 

5.1. Model parameters 



Table 3. Summary of the recovered model parameters. Chi- 
squared values refer to the "best-fit" model. Reported values are 
the median and the 16% and 84% quantiles around it from the 
marginalised one-dimensional posterior distributions. Results 
from the two priors cases are Usted for comparison. Derived 
constraints on additional secondary parameters are also listed 
(slopes are defined in Sect. 5.2). 
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4.3. inclusion of more informative priors 



We consider two working hypotheses for the modelling of 
SL2SJ02176-0513. As our primary and default strategy, we 
adopt the most informative prior - based on the stellar mass de- 
rived from the u*g'r'i'z' aperture photometry - i.e. a uniform 
prior 2.1 < M. < 4.1 x lO^M© (hereafter refered to as M, 
prior). As a sanity check, and to investigate the effects of the 
informative prior, we also consider the broad uniform priors 
on the primary parameters <Tg,p, edm, PA^m, r^Am, CTdm and M, 
introduced above. This case will be refered to as broad prior. 
Although we will refer primarily to the parameters inferred from 
the Mt prior, we will also report for completeness the parame- 
ters inferred from the broad prior. As a further check, we inves- 
tigated the assumption that the dark matter halo and the stars 
are aUgned on the scales probed by the strong lensing features 
(6 < /? < 12 kpc). It turned out that this prior provides virtually 
the same constraints as the broad prior case, and therefore we 
will not discuss it, for the sake of conciseness. 



We show the marginalised posterior distributions for both 
the M, and broad prior cases in Fig. 4. We illustrate pair-wise 
degeneracies between the most relevant parameters by plotting 
the complete (except for the unconstrained core radius dimen- 
sion) set of 2D marginalised posterior PDFs. One-dimensional 
marginalised posterior distributions are also shown in the diag- 
onal panels with overlaid, in blue, the implicit prior that we as- 
sumed for the derived parameters when making the broad prior 
assumptions about the primary parameters. This shows that - by 
construction - our apparently uninformative broad prior model 
assumptions will tend to maximise the mass component that is 
in the form of stars. 

We note the following about our inferences: 

- The misaUgimient between mass and Ught APA = 2Q^^^ deg 
is preferred by the M, prior. However, the same prior also 
prefers a rounder halo, making the misalignment of little 
physical significance. With the broad prior, the misalignment 
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Fig. 4. Confidence regions for all model parameters, except the lens galaxy core radius, which is unconstrained. For the off-diagonal 
2D plots, contours represent 1 and 2-cr confidence levels. Along the diagonal, ID marginalised distributions are shown. Two 
different prior PDF assignments are considered, the basic broad prior (black), and the prior (red) based on the {u*g'r'i'z'} -hand 
CFHTLS photometry and for which stellai" mass M, = 2. 1 - 4. 1 x 10' 'Mq. The dot-dashed blue curve in the /dm(< ^e) PDF panel 
represents the (un-normalised) implicit broad prior on this derived parameter 



of the dark matter halo and the stellar component looses its 
statistical significance as well. 
- The lens model requires a substantial contribution from the 
group halo: for the M» prior we infer (Tgjp = 550;^24okms"' 
which is consistent with both the velocity dispersion of 
galaxies in the group (cr^^'^ = 770+ 170kms"') and the pre- 
diction from the X-ray luminosity scaling relation {cr^^ - 

520 + 120 km s"'). We observe that the more massive the 
group, the less needed is the contribution of a dark mat- 
ter halo on the scales probed by the Einstein ring. A less 
massive group halo is more in line with M, prior as we get 



cTgip = 690^^45 km s~' for the broad prior The informative 
prior thus allows us to obtain tighter constraints on the mass 
of the group. 

- The overall group contribution to the mass budget within the 
Einstein radius is about 4%, while we find it to contribute 
about 8 - 13% to the local density at the Einstein radius. 
The implications for the mass budget in the lensing galaxy is 
thus not negligible. This important issue is discussed in more 
detail below. 

- The galaxy halo dark matter fraction inside the effective ra- 
dius is highly sensitive to the choice of prior For the M, prior 
we obtain fdm(< Re) = O-^l^^oog- Assuming the broad prior 
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would give a much lower value of /dm(< ^e) - 0.06^^ ^g. The 
number obtained with the informative M» prior is consistent 
with local estimates from the SLACS survey (e.g. Gavazzi 
et al. 2007; Bolton et al. 2008b), where typically ~ 30% of 
the projected mass inside the effective radius is found to be 
in the form of dark matter, independent of any assumption on 
the stellar populations. The broad prior yields a much smaller 
number, which is hard to reconcile with other observations. 
This lends confidence to the use of the informative M, prior 
as our primary choice. 
- The core radius in the dark matter distribution is uncon- 
strained by the data and therefore we do not show the pos- 
terior distribution, which is indistinguishable from the input 
prior. 

5.2. The mass density profile of SL2SJ021 76-0513 

Thanks to the appearance in SL2SJ02 176-05 13 of two sets of 
multiple-images from two sources at different redshifts, we can 
expect this lens, by comparison with the recently discovered 
double Einstein ring SDSSJ0946+1006 (Gavazzi et al. 2008), 
to provide tight constraints on the inner parts of its total den- 
sity profile. Gavazzi et al. (2008) showed how important it is to 
take into account the multiple deflector plane ("compound lens") 
nature of such systems, since the nearest background source fur- 
ther bends the light rays coming from the most distant source. In 
the case of SL2SJ02176-0513, the two sources are much further 
away from the foreground galaxy than in SDSSJ0946+1006, 
where the inner ring is at z = 0.609. Furthermore, the source 
positions are not as closely aligned with the lens galaxy as in 
SDSSJ0946+1006, where the sources gave rise to almost com- 
plete Einstein Rings. Thus, the impact parameter of the most 
distant source in the plane of the intermediate source is much 
larger than for SDSSJ0946H-1006 in units of the Einstein radius 
of the intermediate galaxy (assuming a typical size galaxy). It 
is therefore not necessary to treat SL2SJ02 176-05 13 as a com- 
pound lens. In addition, the image configuration provides too lit- 
tle information to put constraints on the mass of the intermediate 
Zarc = 1.847 galaxy, confirming that its effect is insignificant for 
this problem. 

The nearby galaxy group will provide a substantial amount 
of convergence at the lens location, making the issue of the pro- 
jected density profile more complicated than for isolated lenses. 
In the previous section we found that about 4% of the mass 
enclosed in the Einstein radius is supplied by the group halo. 
Locally, at the Einstein radius, the group halo accounts for ~ 8% 
of the surface mass density. Hence the results one gets for the 
density profile slope of the lens galaxy strongly depend on the 
external convergence. This issue is well illustrated by Fig. 5, 
which shows the recovered surface mass density profile: at each 
radius the density is evaluated from the median density over all 
MCMC samples (the envelope shows the 68% confidence inter- 
val, estimated with the 16 and 84 percentiles). In black is shown 
the total density profile - including the stars, the DM halo and 
the (approximately) constant convergence from the group. 

The red curve shows the contribution of stars, seen to domi- 
nate the density profile on small scales as expected. 

Let us now quantify the density profile slope more precisely, 
focusing on the logarithmic projected surface mass density slope 
evaluated at the effective radius: 




radius [kpc] 

Fig. 5. Total surface mass density profile. The effective radius 
and the Einstein radii are shown as vertical dashed blue and dot- 
dashed green lines (respectively at ~ 6.2 and ~ 10.3 kpc). The 
envelope around median values shows the 68% confidence in- 
terval. We show the median stellar density profile (red curve), 
the median DM halo profile (green) and the median eff'ective ex- 
ternal convergence (blue horizontal line) that successively dom- 
inate from left to right. 

This slope parameter would be - 1 for an isothermal density pro- 
file. We can compute a for every MCMC sample and thus build 
up the posterior PDF P(Q'|data). We plot this distribution in the 
top panel of Fig. 6; again, we are showing here the total profile 
slope, including dark halo, stellar and group mass components. 

Referring to Fig. 5, we can see that out to ~ 4 eff'ective 
radii the total density profile (black) is well-approximated by 
a single power-law. Fig. 6 shows that this power law index is 
a = -1.03^g jg, consistent with the isothermal case. However, 
the bottom panel shows the posterior probability distribution for 
the slope parameter a' obtained when we remove the external 
convergence from the mass budget. This shifts the result towards 
a slightly steeper density profile a' - -1.1 l^JJ j^. Constraints on 
a and a' are listed in table 3 for the broad and M» priors. 

As in the case of the dark matter fraction /dm(< ^e), the 
constraints we infer on these derived (secondary) parameters 
need to be compared to the input priors we assumed. The blue 
curve in Fig. 6 shows the un-normalised implicit prior proba- 
bility distributions for a. This turns out to be fairly uniform, 
with only a broad peak close to -1 arising from our use of the 
pseudo-isothermal model DM halo (the high Sersic index stellar 
mass profile favours slightly steeper-than-isothermal total den- 
sity slopes at Rg). 

The difference between the results obtained with the two pri- 
ors can be summarised as: 

- the steep density slopes are mostly driven by the high pre- 
dicted stellar mass with the broad prior; 

- a lower stellar mass, as favoured by the informative M, prior, 
implies density slopes that are closer to isothermal with a - 
-1.03!»:^2anda' = -l.ll!»:2^; 

- the slopes a and a' get closer to one another, reflecting the 
reduction in the contribution to the mass budget from the 
nearby group when the stellar mass is lowered (and the more 
extended DM halo is increased in mass to compensate). 
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Fig. 6. Probability distribution for the logarithmic surface den- 
sity slope a at the effective radius location for the total (stel- 
lar+DM+group convergence) density profile. The colour cod- 
ing for the curves showing posterior probability distribution with 
different priors is the same as in Fig. 4. 

5.3. The external convergence - mass profile degeneracy 

Even with two multiply-imaged systems at two different red- 
shifts as we have here, gravitational lensing alone cannot tell us 
both the density profile slope and how much mass is in the form 
of a uniform external mass density. However, the degeneracy is 
broken here, via the constraints we have placed on the associ- 
ated external shear provided by the group halo, the assumption 
that the group halo has an isothermal mass distribution, and the 
functional form assumed for the mass density profile of the lens 
galaxy (because the two multiply-imaged systems are not ex- 
actly at the same projected location in the sky). For example, by 
making these assumptions, the external convergence and shear 
strength are fixed to be equal. 

As we show in Appendix A, one could also break this degen- 
eracy with a measurement of the stellar velocity dispersion of 
the lensing galaxy; this would allow us to infer both the density 
slope and external convergence at once. We postpone this analy- 
sis for future work. We note that the combination of stellar kine- 
matics of the lens galaxy in a single aperture and strong lensing 
constraints coming from a single multiply-imaged system {e.g. 
Miralda-Escude 1995; Sand et al. 2004; Treu & Koopmans 2004; 
Koopmans 2006; Koopmans et al. 2006; Czoske et al. 2008) can- 
not break the degeneracy between both density slope and exter- 
nal convergence by itself: additional information is necessary. 
This could include for example spatially resolved stellar kine- 
matics, information on the external mass distribution, or mul- 
tiple lens planes together with assumptions on the form of the 
mass distribution as in the case presented here. 

6. Discussion & Conclusion 

In this paper we have attempted to produce a detailed lens mod- 
elling of one of the most complete systems in the growing 
"SL2S" sample of strong lensing events in the CFHT Legacy 
Survey. SL2SJ02 176-05 13 is a partial Einstein ring in the out- 
skirts of an X-ray-emitting group of galaxies. The observational 
dataset available for SL2SJ02 176-05 13 is quite remarkable, con- 
sisting of deep u*g'r'i'z' CFHTLS photometry, high resolution 
data from HST/ACS snapshot imaging, a published XMM X-ray 



map, the galaxy velocity dispersion of the host group (based on 
18 radial velocities; Geach et al. 2007), and a spectroscopic mea- 
surement of the deflector and source redshifts from Keck/LRIS. 

In addition, the system exhibits an extra double-image sys- 
tem, which we have determined to lie at higher redshift than 
the initially detected arc/counter-arc system. This system gives 
an important additional handle on the radial mass distribution. 
To our knowledge, SDSSJ0946-H1006 (Gavazzi et al. 2008) and 
SL2SJ02 176-05 13 constitute the only two known galaxy-scale 
systems with multiply-imaged sources at different redshifts. 

The influence of the nearby group on the lens configuration 
is detected as an external shear that accounts for the slightly mis- 
aligned overall potential and stellar component. The lens model 
including a prior on the stellar mass from the CFHTLS pho- 
tometry seems to favour a massive group with velocity disper- 
sion (Tgfp - 550;^ km s"', consistent with X-ray luminosity- 
predicted (Geach et al. 2007). The group environment contribu- 
tion is somewhat degenerate with the mass of the dark matter 
halo of the lens galaxy: the higher the group mass, the lower the 
dark matter fraction in the deflector galaxy. Similarly, if the shear 
supplied by the group is reduced, the galaxy dark matter halo and 
the stellar component must be preferentially misaligned. 

Despite these degeneracies between model parameters, the 
total radial projected mass profile is well-constrained by the dou- 
ble source plane nature of SL2SJ02 176-05 13. However, we have 
showed that the slope of the galaxy stellar-plus-dark matter mass 
density profile at the effective radius differs according to whether 
the convergence supplied by the group is taken into account in 
the mass budget or not. The effect the group has on the apparent 
density profile is to make it shallower by 8% (from a' - -1.11 
without the group contribution to a - - 1 .03 with the group con- 
tribution). Again, this relatively steeper than isothermal profile 
is consistent with the lens PGl 1 15+080 in the environment of a 
group (Treu & Koopmans 2002). Statistical analyses of SLACS 
lenses also show that they tend to have slightly steeper density 
profiles when found in denser environments (Auger 2008) or, 
more precisely, when they are satellite of massive halos (Treu 
et al. 2008). This latter result is also supported by simulations of 
Dobke et al. (2007) and Limousin et al. (2007b). 

As mentioned in the Appendix, a velocity dispersion mea- 
surement of the stars in the central galaxy would help to break 
the remaining degeneracy in this system, tightening the con- 
straints on the galaxy density profile and the group mass. We 
are currently pursuing this observational effort, and we are now 
working on the spectroscopic determination of the velocity dis- 
persion of stars in the lensing galaxy. Likewise, the continuing 
growth of the SL2S sample will allow us to tackle the problem 
of understanding massive galaxy mass distributions with better 
statistics. 
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Appendix A: Constraints on density slope and 
external convergence with three aperture 
masses 

In this appendix, we demonstrate that the combination of two 
multiply-imaged sources at different redshifts, which define two 
projected aperture mass measurements, and an additional mea- 
surement of a small-scale stellar velocity dispersion, that can be 
approximated as a three-dimensional aperture mass estimate, al- 
low us to constrain simultaneously the density profile slope pa- 
rameter and the external convergence due to mass lying both 
along the line of sight and in the lens plane outside the lens 
galaxy (but that does not affect the lens galaxy's stellar dynam- 
ics). 

Let us define a density profile of the form: 



-r 



P(r) = Poir/ro) 

3-r 



(A.l) 
(A.2) 



with 1 < y < 3. Once projected along the Une of sight, we get 
the projected surface mass density: 

m) = 2:o(/^/ro)'-^=PoroV^^^^^^p^-!^(/^/ro)'-^ (A.3) 

r[r/2] 



M2diR) = Sor2(/^/ro)^-^ 

3-r 



(A.4) 



We now define a possible constant surface mass density 2sheet> 
and three mass measurements M^, M^^ and M^j at radii Ri, R2 
and /?3 that satisfy: 



Ml, = poA(y)(Ri/ro)^-^ + nR%t^,t 
Mid = PGA{y){R2ln,f-y + nR%beet 
Mm = poB{y){Riln,f-y . 



(A.5) 
(A.6) 
(A.7) 



We can readily see that only in the case of three observables 
can one hope to constrain all three of po, y and Ssheet (''o be- 
ing an irrelevant rescaling here). In the equations above we 
have substituted A(j) = rl2K^^^r[(y - l)/2]/r[7/2]/(3 - 7) and 
B{y) = rg4;r/(3 - y) for brevity. Schematically we thus obtain: 



r[(y-l)/2] 

r[r/2] 



2 R'i'RlMii- 



■R\Mli 



■^sheet — 



^|^MM RlR\-y - R\r\-^ ' 

im\,r\-'-m1^r\-' 



■RiR'r 



(A.8) 



(A.9) 



where Eq. (A.8) allows us to solve for y and then Eq. (A.9) yields 

^sheet- 



